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We predict an existence of an new dissipation channel leading to attenuation of ultrasound in 
disordered conductors and superconductors with perfect electroneutrality. It is due to slow diffusion 
of thermal energy. We show that in doped silicon ultrasound attenuation may be enhanced by 
a factor about 100. Similar effect is also studied for s-wave and d-wave superconductors. The 
latter case is applied to BSCCO family where strong enhancement of ultrasound attenuation is 
predicted. For usual s-wave superconductors new dissipation channel might be important for very 
low-electron-density materials near the BCS-BEC crossover. 

PACS numbers: 71.10.-w, 62.80.+f, 74.25.Ld 


I. INTRODUCTION AND MODEL. 

Ultrasonic attenuation in metals have already been 
studied for a long time^^® and may seem to be fully un¬ 
derstood. The ratio of ultrasound attenuation rate aiw) 
to the sound frequency uj in clean metal is small due to 
adiabatic parameter mkp/pm 1 where m is the elec¬ 
tron mass, kp is the Fermi wave vector and pm is the 
density of the material. The simplest model for electron- 
phonon interaction is the scalar vertex Frohlich model 
with an extension due to Migdal®’^. This model is ap¬ 
plicable for clean metals. However, it is well known that 
the Frohlich model is not adequate when phonon wave¬ 
length 27r/5 exceeds elastic electron mean free path 1. In 
this dirty limit, when ql 1, the conventional theory of 
electron-phonon interaction in a disordered conductors 
leads to the Pippard ineffectiveness condition(PIC) that 
tells that the ultrasonic attenuation at small wave vectors 
is suppressed by a factor ql 1^. 

The PIC phenomenon results from strong Coulomb in¬ 
teraction that prohibits any charge disbalance in the sys¬ 
tem (perfect screening condition). Still some mechanisms 
were found recently which increase inelastic electron- 
phonon rate even in presence of perfect screening. The 
first of the them is present in a multiband electronic 
spectrum®’®, the second one is realized when impurities 
do not quite follow a motion of a lattice^®. Still another 
possibility is related with a deviations from the perfect 
screening condition that become prominent in a strongly 
disordered conductors, where kpl is not very large®’^^. 

In the present paper we show that even in the sim¬ 
plest case of one-band electron spectrum and a perfect 
screening, the conventional PIC theory might still be un- 
sufficient. The reason is that electron liquid possesses 
an intrinsic diffusion mode that is present even under 
the condition of strict electroneutrality. It is an energy 
diffusion mode that plays a role similar to the spin po¬ 
larization mode studied in Ref.®. Coupling of phonons to 
this diffusion mode results in an additional attenuation 
of longitudinal sound waves. 




Figure 1: (Color online) Phonon self energy for (a) local pro¬ 
cesses and (b) diffusive channel, when phonon converts into 
diffusive mode. Throughout this paper we use Keldysh dia¬ 
grammatic technique. 

We consider a dirty conductor with an electron action 

Sei,n = J dt{dr)tjj* [idt - (^(p) -f U{r))] ifi (I) 

with U being disorder potential and i =t, i correspond¬ 
ing to spin indices. Electron-phonon interaction is given 
by 

Sei-ph,n = J dt{dr)tp* [{ppVp/d)daUa{r)- (2) 
-aa(u„(r){7(r))] 'ifi, 

where the first term describes interaction due to the mod¬ 
ulations of ionic density while the second term corre¬ 
sponds to distortions of disorder potential, the latter be¬ 
ing described by a correlator 

{U{v)U{y')) = u5{r — v'). (3) 

It can be shown, using comoving frame of 
reference®^®’®’^"^, that the action in the form of Eq.(2) is 
equivalent to 

Sei-ph,n ^ J dt{dv)ip*T‘ff {dpUafifi, (4) 
with the vertex 

T’ff = [paVp - {ppvp/d)5c,ii] , (5) 

where pp, vp are Fermi momentum and velocity, d is the 
dimensionality of the electron system. 




II. ULTRASOUND ATTENUATION IN 
NORMAL CONDUCTORS. 
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Local processes Attenuation rate a(uj) is related with 
the imaginary part of the phonon self energy, 


a(cLi) 


^ ImS'^(w, q) 
Pm^ 


( 6 ) 


Local contribution to ultrasonic attenuation is due to 
processes which involve creation and annihilation of elec¬ 
tron and phonon states on the spatial scale comparable 
to the elastic mean free path. Such a process is depicted 
by the diagram in Fig.2 a) and leads to a well known^^^ 
result (it can be obtained via the calculation of the dia¬ 
gram Fig. fa with the vertices defined in Eq. 5): 


anj=2ci'^^Dq^ (7) 

Pm 

where u is electron density of states per one spin, D is 
diffusion coefficient and ci = 2{d — l)/d{d + 2) is just a 
numerical coefficient. 

Energy diffusion processes. For local processes infor¬ 
mation about direction of electron motion is retained be¬ 
tween the acts of absorption and emission of phonons, 
thus the averaging of the product of two vertices over di¬ 
rections of momentum in the diagram (Fig. fa) leads to 
a large result: 

^r„(p)r„(p)^^(8) 

However, when diffusion processes are allowed, the 
same information is lost between absorption and emis¬ 
sion of phonons. On a formal level, it is seen from the 
diagram (Fig. lb) where an impurity ladder is inserted 
into the diagram disconnecting the electron-phonon ver¬ 
tices. Multiple collisions with impurities hold electrons in 
the regions of phonon absorption and emission leading to 
independent averaging of the interaction vertices F over 
directions of momenta. That results in a much smaller 
values of effective vertex, which is now of the order of the 
temperature T <C pf'Cf, 



Figure 2: (Color online) When coupled to diffusive mode, 
electron vertex is always averaged over the momentum p that 
goes in the block of two Green function preceding the diffuson. 


where subscripts ± stand for (e ± a;/2, p ± q/2) respec¬ 
tively. In general, effective vertex depends on both en¬ 
ergies e± = e ± w/2 but in the limit uj e ^ T) depen¬ 
dence on u! is negligible. New vertex Tph-d,ni£+, £-) is 
then given by 


(F)f = (11) 


where 


d{pv/d) _ ( pFVFd\n.v\ 

Diffusive electron modes evolve on a much larger 
timescales ~ {Dq^)~^ implying a more effective dissi¬ 
pation, despite weak phonon-diffuson conversion. While 
the square of the vertex is smaller by the factor (T/etf)^, 
the dissipation due to diffusion processes is enhanced by 
a factor {ql)~^ (at the lowest frequencies enhancement is 
saturated by a factor v\ls^, s being a sound velocity). 

To the best of our knowledge, such processes for a 
single-band conductor under the condition of the perfect 
Coulomb screening had always been neglected so far. The 
corresponding diagram is shown in Fig. lb. Its evalua¬ 
tion leads to the following result (see also Eq.(6)), valid 
aX Huj T: 






dfie,T) \ 
de ) 


vO 


w), 


(13) 


(r„(p))^^e~T (9) 

This is the reason why diffusion modes are usually ne¬ 
glected in electron-phonon interaction in disordered con¬ 
ductors. 

It is convenient to define effective electron-diffuson ver¬ 
tex that depends on electron energy instead of a quickly 
relaxing momentum. Diagrammatically (see Fig. lb) 
there is always a block of two Green functions that stands 
in between the phonon and diffusive mode and results in 
averaging over the electron momentum. Therefore an ef¬ 
fective phonon-diffuson vertex (F)^ may be defined with 
the help of the following integral representation: 

(r)f (^) = ^/[r“^(p-,p+)] Gt (10) 


where Q:|,^^(a;) is the partial contribution of electrons with 
energies in the interval (e, e -|- de), 

^ {2vReV{uj,q)) . (14) 

Pm 

Eq.(14) contains diffusion propagator V{uj,q) (the corre¬ 
sponding diagrams are shown in Fig. 3) equal to 


V{uj,q) 


1 

—iuj + Dq'^ ’ 


(15) 


where r = Z/uf is an electron elastic scattering time and 
the diffusion coefficient D = rvp/d. (Note that r = 
(2Trh'u)~^ with electron density of states f and u defined 
in Eq.3) Calculation of the integral in Eq.(13) leads to 
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the following result for the ultrasound attenuation rate 
at frequency uj and temperature T: 

27r2 vD 

^ Pms^ 1 + {Doj/s^)^ ■ 

It is useful to present it also in the form of the ratio of 
(16) to the ’’local” result, Eq. (7): 


an,d _ ^ f . 

anj 3ci \Pfs ) + a;2 ’ 



The two mechanisms of dissipation lead also to dif¬ 
ferent dependencies of attenuation on temperature and 
frequency: ’’local” attenuation rate is only weakly T- 
dependent but grows as On the other hand, attenua¬ 
tion due to ’’diffusive” mechanism is proportional to T^. 
Dissipation due to diffusive mechanism is frequency- 
indepedent at high lo tOc = jD and goes as at 
low frequencies. Formulae (16,17) can be used in order 
to extract the value of the electron diffusion constant D 
from the data on the ultrasound attenuation. 

For thermal phonons with hw ^ T attenuation due 
to diffusive channel is always small, an,d ctn,i- How¬ 
ever, for the ultrasonic attenuation at low frequencies 
u! ^ T/hihe situation can be quite different, especially in 
doped semiconductors at moderate temperatures. Con¬ 
sider, for example, heavily doped Si with n = 10^°cm“^, 
m = 0.36mo, ppl — 10, s « 8 • lO^cm/s. At the temper¬ 
ature T — O.lEp ~ 200Ar one finds 


« 100, / = 27rw < lOGHz (18) 

^n,l 


i.e. attenuation is enhanced by two orders of magnitude 
due to contribution of the diffusive channel. 


III. S-WAVE SUPERCONDUCTORS 
A. Formulation of the model 


We consider BCS superconductor (SC) with an s-wave 
pairing and microscopic interaction constant g, 


Sei,sc-int= J dt{dr)[gtjj*{r)i;*{r)-tfj{r)i;{r)]. (19) 

In the case of a dirty superconductor it is convenient 
to define Green function as 

Gf^(e,p)=-*(M'„^^) (20) 

= [(£ -f i0)f3 - ^fo - Aiif2)]~p, 

where 4-component spinor 4^ is defined as follows: 


■Kp) = ^ 


f V’pt \ 


= 5-+f3. (21) 



Figure 3: (Color online) Diffusive modes correspond to im¬ 
purity ladders. The quantity !?(...) in the text is given by 
the sum of this ladder, excluding one overall factor 

coming from the first impurity line. 


The 4-dimensional Nambu-Gor’kov space corresponds to 
the product of spin and particle-hole spaces. However, 
for our purposes only particle-hole subspace is relevant, 
as is evident from Eq.(20), where only particle-hole Pauli 
matrices f are present. We will thus ignore the spin struc¬ 
ture, effectively working with 2-component fk, 'k spinors. 

For a disordered s-wave superconductor self-consistent 
Born approximation gives^® 


G^{£, p) = V{s) (sfs - A(*f2)) - $fo 


with 


( 22 ) 


77 (e) = 1 + 


2tE 


E = - AA (23) 


B. Energy density mode 

1. Diffuson 


Diffusion modes are described by impurity ladder 
shown on the Fig. (3). This ladder can be analyzed with 
the help of a Bethe-Salpeter equation 

{'D)-\e,uj,q) =l-E{e,uj,q) (24) 

with the self-energy 

E{e,oj,q)=u J (dp) G^{e-,p) ig: G^{e+,p). (25) 


Green function is a 2 x 2 matrix in the particle-hole space 
(2-dimensional spin space is irrelevant, see Eq.(21)). 
Thus, V and S are matrices in a 4-dimensional space 
constructed as a product of two 2-dimensional particle- 
hole subspaces. These 4x4 matrices D and 2 could be 
interpreted as superoperators acting on 2 x 2 operators 
such as effective electron-phonon vertex T (for an s-wave 
state r is defined further in this chapter, see Eq.(30)). 
For example, if A, B and X are 2x2 matrices, then 
Y = AXB is as well a matrix of this type. This means 
that A® B describes a mapping A —)• F thus indeed 
being a superoperator. 

At zero external frequency and momentum the self¬ 
energy reads as follows: 


^ 1 


{efs - A ( 77 - 2 )) G (£f3 - A ( 7 - 72 )) 

2(e2 _ ^2) 


(26) 
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Analysing this matrix structure one finds two massless 
modes corresponding to operators 

To and £:f 3 -A(*f 2 ). (27) 

The To mode is symmetric in particle-hole space and 
corresponds to charge density, while the asymmetric 
mode £73 -I- A ( 172 ) corresponds to energy density. The 
corresponding term in the matrix propagator defined in 
Eq.(24) is given by 


Here X = vg is a dimensionless BCS coupling constant. 
Below in this section we consider temperatures and all 
other relevant energy scales to be well below the gap 
T, w, < A. 

Variations of A arise either directly through the density 
of states or through interaction constant g, 

dlnA dlni' ding 
dlnp dinp dinp 


{V){e,u},g) 


2'kvt'^ 




where 




1 

[-i(A+ -E_) + Dq^~] 


(28) Within our model dlniy/dinp = {ppvp/d){dlni^/dE)] it 
arises due to the shift of the chemical potential in the 
presence of acoustic wave. 

Defining effective vertex (A) analogously to Eq.(30) 

(29) 

we get 


and E± = \J{eE oj/2Y — A^. Charge density fluctua¬ 
tions are prohibited by strong Coulomb interaction, so 
we should not consider them here. 

Strictly speaking, there are two more diffusion modes 
of the ’’Cooperon” type which are related to the 
charge conversion between condensate and thermal ex¬ 
citations. Contrary to Eq.(25), these modes are con¬ 
structed with matrix products of the type of G^(£_, p)® 
G'^(—£+,p). The corresponding propagator is propor¬ 
tional to (—i(i?+ -I- E_) + Dq^)~^. It is important that 
in the limit Dq^ -\/TA the contribution from these 
modes is much smaller than that of the energy diffusion 
channel, so we neglect Cooperon channel in the following. 

2. Effective vertex 


(A), (£) = >fA (^^3 - A(if2))<5„/J. (35) 

The effective vertex (A)^ is substantially enhanced due 
to singular density of states, in contrast with the vertex 

(f).- 


C. Ultrasound attenuation due to energy diffusion 

Similarly to the normal metal case, contribution of dif¬ 
fusion channel to ultrasonic attenuation is given by the 
following integral (assuming uj T): 


When defining equivalent phonon-diffuson vertex we 
now have to pay attention to the matrix structure in the 
Nambu-Gor’kov space (f„ = r„ • 1), 

(f)j£)= y'(dp)G«r„(p_,p+)G^, (30) 

that modifies the vertex (11) into 


asA^) = I ai%) (36) 

0 

where is the partial contribution of quasiparticles 

with energies in the range (£, £ -I- ds), 


(f)^ (£) = xr(^£f3 - A(7f2))(5a;3. (31) 


We see that there is indeed no coupling between ultra¬ 
sound and charge density mode. 

Generally, acoustic wave modifies effective BCS cou¬ 
pling constant A leading to an additional electron-phonon 
vertex. This vertex has the same matrix structure as the 
order parameter: 


A = XrA(Afi)(5a/3 

with constant >«’a being equal to 

dA 1 /din A A 

^ dlnp BCS A ydlnpy 


(32) 


(33) 


^ 9)) > 

(37) 

(w, q) being a diffuson in a superconducting state 
given by Eq.(29) and E = y/e'^ — We have also added 
subscript (n) to normal metal DOS Vn to avoid possible 
confusion. Note that the contribution due to variations of 
the effective BCS interaction constant (described by the 
vertex (A)) is enhanced due to singularity in the density 
of states at the gap edge, iy{e)/i'n — A/^/e'^ — A^. 

Substitution of Eq.(37) into Eq.(36) and integration 
yields final result for ultrasonic attenuation due to energy 
diffusion channel in s-wave superconductors: 
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Q!s,d(w) = 2- - exp 

PmD 


A 

T 


f (A/r) In + 4>!rxrAAr + 4>ir^T2^ 

1^2 (xr^A^ + 4xrxrAAT + Sxr^T^) x ^ 


W ^ Wc 
W <C Wc 



(38) 


where A(a;) = w + A(wc/w)^ and Wc = s^/D is the 
crossover frequency in the normal state. Equations (38) 
were derived under the conditions 


w,£>g^<T<A<r ^<^pfVf (39) 


and contain two characteristic crossover frequencies, 
Wc.si = a;c(A/T)i/2 , and uJc,s 2 = Wc^^A^/^. 

Note that ujc,s 2 describes a weak logarithmic crossover 
that results from the DOS singularity in the supercon¬ 
ducting state; (A/T) factor in the first line of Eq.(38) 
is the result of DOS behaviour as well. This crossover 
exists only if uJc,s 2 > ^c,si or equivalently only for tem¬ 
peratures T > (Awg)^/^. The frequency of the second 
crossover uis^c 2 can be small or large in comparison with 
Ws,cij depending on specific material. 

To compare attenuation due to the energy diffusion 
channel (38) with the one produced by usual local pro¬ 
cesses in superconducting state (we denote it as a^,;) note 
that the latter is proportional to the density of normal 
electron-hole excitations^®. Thus the ratio 


as,i 

C^n,l 



dm 

de 


2 exp 


A 

T 


(40) 


so that ag^i oc as in the normal state. At high fre¬ 
quencies uj ^ uJc,s 2 we thus find that the role of energy 
diffusion channel grows with w decrease: 


Qs.d 

as.i 


i_!_/ 2^1 T 

2ci \ ^ T A(a;) 




+ 4.>AT)x(^)^ 


(41) 


In this frequency range attenuation due to energy dif¬ 
fusion has only weak logarithmic frequency dependence; 
it also depends on temperature in nontrivial way that 
does not reduce to the T-dependence of normal electron 
density. 

At lowest frequencies, when w < Wc,si, frequency de¬ 
pendence of tts is of standard form, but temperature de¬ 
pendence differs from that of local contribution: 


pairs”. Particular example of this kind is presented 
by recently discovered heavy-metal compound YPtBi^^ 
with conduction electron density at temperatures about 
Kelvin range as low as n = 2 • With the values 

of other parameters taken from Ref.^^ as m = O.IStoq, 
electron mean free path I = 130nm, = 0.77K, lof) = 

200Ar and sound velocity^® s = 2 • 10®cm/s, we find for 
T = 0.2K the ratio of the energy diffusion contribution 
to the standard PIC result as,i'. 

^ « 1, f = 2 ttu< IGHz. (43) 

IV. D-WAVE SUPERCONDUCTORS 
A. Hamiltonian 

We consider here strongly a model of anisotropic d- 
wave superconductor and neglect electron dispersion in 
the direction transverse to the layers. Then the order 
parameter depends on direction of momentum(Fig.4) in 
the following way: 

A(p) = Ao(cospa;a — cospytt), (44) 

where a is the lattice constant. At low temperatures 
T «C Ag this leads to angular dependence of ultrasonic 
attenuation due to local processes. At the same time, the 
contribution of diffusion channel is isotropic (see Fig. (5)). 

For d-wave superconductor it is convenient to define 
Green functions somewhat differently from s-wave case, 
namely: 

= [(£ -k )i 0 )f 3 - ^ - (ff 2 )A]-\ 

145) 

see the difference with Eq.(22), where we had used = 
f 3 'I'+ instead of 4'+. Energy of excitations above the 
d-wave SC state vanishes at 4 nodal directions in the 
momentum space of the order parameter, see Fig. (5). 
Near these directions the excitation spectrum can be lin¬ 
earized: 


+ 2hx^AT+ 8xm). (42) 

as,i Cl pm ^ ^ 

Inspection of the ratio (42) shows that usually it is 
rather small; interesting exception is presented by super¬ 
conductors with extremely low electron density, which 
are not far from the crossover to the regime of ’’local 


m VFki, A = Vgh, (46) 

where momentum k is measured from the node, k = 
P Pnode, the basis is individual for each node (see the 
Fig. 4). Due to the presence of low-lying excitations, all 
thermal effects at T ^ are described by some power 
laws of T, instead of an exponential behavior in s-wave 
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Figure 4: (Color online) Schematic Fermi surface of a d-wave 
superconductor. Background color represents order param¬ 
eter A(p) angular dependence revealing four nodes on the 
Fermi surface. Each of them has its own electron Fermi 
vi = d^/dp and gap V 2 = dA/dp velocities. We thus choose 
four different coordinate systems (fci, fc 2 )*'*Hfor clarity only 
one is shown) around each node with axes aligned with local 
electron Fermi and gap velocities, x, y directions correspond 
to [ 100 ], [ 010 ] crystalline axes; a is the lattice constant. 


(in units of jh). Self-consistent Born approximation is 
valid as long as iC ^ 1 ; since the ratio vp/vg is a large 
parameter (according to Ref.^°, it is about 20 for BSCCO 
superconductors), thus normal-state conductance should 
be sufficiently large, Gn > 100 , for this approach to be 
valid. 

Asymptotic behaviour of 7 (e) and Z^e) is of particular 
interest, it is given by 

£ < 70 : 7(0) = Aoe-i/^, Z(0) = K, (53) 

(«) 

where 70 = 7(0) is a scattering rate at the Fermi surface. 
Similar results were obtained previously for graphene^^. 
Below we restrain ourselves to the intermediate temper¬ 
ature range 70 T <C Aq. 

To describe electron-phonon interaction, it is still 
sufficient^"^ to start from the electron stress tensor Eq.(5), 
but now a momentum-dependence of the order parameter 
A(p) should be taken into account: 

f = PadpH - (55) 

^ /'{vfPf/ 2)+VFkF (pf + kF)vg{if2) \ 

\ kgVp -{pFVF/2)+kgVg{if2)) 

where we have used H = + A(ff 2 ). 



case. Now we should take into account random disorder 
potential, which we consider in the white-noise limit: 


B. Ultrasound attenuation due to local processes. 


{U{r)U{r')) = u6{r — r') (47) 

Scattering by disorder results in renormalization of Green 
functions, leading both to a scattering rate 7 and to 
renormalization of the energy: £—)■£. In the self- 
consistent Born approximation we find 

{G^p{s,P)) = [{e + ij)f3 - ^fo + {if2)A]-^, (48) 

where e is related to a residue of the Green function, 

£ = Z{£)e{e). (49) 


Quantities 7 and £ are to be determined by self- 
consistency equations^® 


I ^ 

1 = K \ — arctan —h In 

\7 7 


-\/£2 -I- 72 


7 7 

where IF is a dimensionless disorder strength, 

ZTTVpVg Vg 


(50) 

(51) 

(52) 


Here iV is a number of valleys {N = 4), and G„ is the 
dimensionless conductance per layer in the normal state 


Like previously, attenuation due to local processes in 
a d-wave superconductor can be depicted by the elec- 



0 n/2 7t 3n/2 2 n 


0, rad 


Figure 5: (Color online) Schematic angular dependence of 
ultrasonic attenuation for local (a;) and diffusive (o-d) pro¬ 
cesses. Q; exhibits strong anisotropy with extremal values 
iai)max showu by horizontal dashed lines. In con¬ 
trast to local processes, attenuation in diffusion channel is 
isotropic. As a reference direction 0 — 0 we choose [110]. 
Minima of ultrasonic attenuation are observed in ±[100] and 
±[010] directions, while maxima in ±[110] and ±[1-10]. The 
angular dependence of the gap A(k) is shown as well (dashed 
blue line). 
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tron bubble diagram (Fig. la), the corresponding analyt¬ 
ical expression reads as 


ad./ = y^Tr [(/+-/_)f(G^ 
where F is an abbreviation for 


G^)f (G^ - G^)] , 
(56) 


f = 




(57) 


with Cq, being a unit vector in the direction of phonon 
momentum qa (while polarization vector of longitudinal 
phonons is essentially the same as 6 ^). In contrast to 
normal and s-wave states, in d-wave state attenuation is 
essentially anisotropic due to two strongly different ve¬ 
locities vpyVg and symmetry of the order parameter. 

Partial contribution of electrons with energies (e, e + 
de) into attenuation rate can be found in the form 


(e) 


:/(^) = 


2TT'^VFVgP„ 


1 -b 


K 




(58) 

where the function F{9) comes from electron-phonon ver¬ 
tices and describes angular dependence of attenuation, 

/ \ 2 


F{0) = sin^ 29 + j cos^ 29. 


(59) 


Integration of Eq.(58) over electron energy (like in 
Eq.(13)) eventually leads to 


ad,i 


F{9) vfpW T 

27r2 Vg PmS"^ 7o ' 


(60) 


Comparing of the result (60) with attenuation in the 
(isotropic) normal state we find 

^ = 27s:F(6l)ln—, (61) 

Oin,l 70 


where K is the dimensionless measure of disorder defined 
in Eq.(52) and 70 is the scattering rate at the Fermi sur¬ 
face, Eq.(53). Note that in the limit T <C 70 the contri¬ 
bution from local processes is very similar to (60), one 
should just replace ln(T/ 7 o) —)■ 2 . 

Attenuation rate Eq.(60) possesses unusual logarith¬ 
mic temperature behaviour oc ln(T/ 7 o) that does not 
correspond to naive expectations for the linear density 
of states, a (X T. The reason is that electron scatter¬ 
ing rate 7 is intimately related to the electron density of 
states; typical scattering rate decreases with T together 
with typical DoS. Namely, Eq.(54) shows that at mod¬ 
erately high temperatures 7 cx T up to a logarithmic 
factor. 


Salpeter equation Fig. (3) 

(5)"^(£,w,g) = 1 - S(e_,£+,q), (62) 

with a self-energy 

S(£,a;,q)=u^ J {dp)G^ G^. (63) 

nodes 

Energy density diffusion mode corresponds to fa eigen¬ 
vector and eigenvalue of the self-energy 

S(£,w,q) j (dp)Tr [faG^^faG;^] . (64) 

nodes 


Summation of the impurity ladder for the energy density 
mode thus gives 

{V){e,uj,c\) ^2u'y{e)V^^\uj,q), (65) 

where is a diffuson propagator in a d-wave state. 


V 




1 


-iuj + 


The diffusion coefficient is 



K + Z{e) 
2-f{e) 




( 66 ) 


(67) 


where (u^) = {vlp + — ^f/ 2 - 

Cooperon modes in a d-wave state are irrelevant at 
temperatures T 3> 7 o- 


2. Effective vertex 

Effective phonon-diffuson vertex is defined in the way 
similar to the s-wave state: 

(f)“''(£+,£_) = u ^ Jidp)G^tf{p)Gt ( 68 ) 

nodes 


Calculation of this vertex requires some caution, see Ap¬ 
pendix for details. Eventually we get, similar to the nor¬ 
mal metal state: 


C. Energy diffusion mode 

1. Diffusion propagator. 

This section proceeds in a way almost identical to s- 
wave case. The diffusion modes are described by a Bethe- 


= e^3^a/3- (69) 

The effective vertex describing coupling to the energy 
density mode turns out to be isotropic oc 5ap once the 
summation over 4 nodes is performed. This happens de¬ 
spite the anisotropy of velocities for each node separately. 








D. Ultrasound attenuation due to energy diffusion 


We obtain the contribution of diffusion processes in a 
complete analogy to the s-wave state: 

aM = ( 70 ) 

where 

^ - • hvd{e)KeV^^\uj,q)\ (71) 

Pm ^ ' 

where Vd{e) = 'y{e)/Tru is a density of states in a d-wave 
superconductor. Then Eq.(71) gives 




7r2 Ug V K ) 1 




(72) 


with a crossover frequency 0Jc,d{£) = We are 

interested only in the case of weak disorder K 1 
and temperatures T ^ jq. Effective crossover frequency 
which enter the result for the total attenuation, as given 
by Eq.(70), is 


7Ts‘^ KT 

^ ^ [1 - ii:in(Ao/T)]3 ’ 7’>7 o- (73) 

In Eq.(73) we have used identity if'ln(T/ 7 o) = 1 — 
7nn(Ao/r). 

At low frequencies the resulting attenuation rate has a 
oc In T temperature dependence 

ad.di^) = 7 — —^ ln(T/7o), w < uJcA^), (74) 

3 Vg PmS^ 

while for higher frequencies the attenuation behaves 
roughly as ocT^, 

Ttt* 

rf(ca) —— 2 T- -j-- tvaTs 5 (7b) 

GOpmvAg [1 - K\n{Ao/T)\^ 

W » UJc,dA)- 

Additional two powers of temperature, oc vs. oc in 
the normal state, are due to nearly linear DOS at high 
energies, J^d(e) oc e up to a logarithmic factors. Note 
that the energy dependence of diffusion coefficient has 
the same origin, oc 7 (e) oc i^d^i^) ^ 

Contribution of energy diffusion channel is especially 
prominent for directions which correspond to minimum 
of the ’’local” attenuation rate (Fig. 5, Eq.60). For such a 
direction at lowest frequencies, when oj < ujc,d, the ratio 
of the contributions from the energy diffusion and from 
local channel is similar to that in the normal case, 


For higher frequencies ui ^ ujc,d for the ratio of attenua¬ 
tion rates we find 



All these results are valid as long as we are in a dirty limit 
ql ^ 1. This condition is equivalent to Dd(T)q^ 7 (T) 
that in turn gives to u 1 q^\t), 

~ ^ ^ tts KT 

° ^ ’~vfAW) 2uF[l-Kln(Ao/r)]5/2- 

(78) 

Using the parameters for BSCCO compounds^^’^^^^"^, 
we estimate the ratio of the (isotropic) contribution of 
the energy diffusion channel (given by Eq.76) to the con¬ 
tribution of local channel at its minimum: 

. « 10, / = 27ra; < IGHz. (79) 

\P^d,l )min 

Here we used the following values of relevant parame¬ 
ters: n = 5 • lO^^CTO”^, vf = 2.5 • lO^cm/s, VFjvg = 20, 
T = lOAT, impurity scattering rate 70 = lAT and sound 
velocity s = 4.6 • 10®cm/s. 


V. CONCLUSIONS 

We have shown that low-frequency phonons in disor¬ 
dered conductors and superconductors experience addi¬ 
tional damping due to coupling between lattice density 
modulations and non-equilibrium distribution of ther¬ 
mally excited quasiparticles. Nonequilibirum distribu¬ 
tion of quasiparticle’s energy then slowly decays lead¬ 
ing to the phenomenon similar to the Mandelstam- 
Leontovich relaxation. The effect is especially strong in 
doped semiconductors at moderately low T/Ef and ul¬ 
trasound frequencies lo T/h. In particular, we have 
estimated ultrasound attenuation in doped Si to be en¬ 
hanced by factor about 100 due to this mechanism. Fre¬ 
quency dependence of the attenuation rate a(w) contains 
typical crossover frequency Wc, see Eq.(17), which de¬ 
pends explicitly on electron diffusion coefficient D. Thus 
measurements of Q:(a;) may be used for determination 
of D. Similar phenomenon exists in superconductors as 
well; for conventional s-wave superconductors it is weak 
usually, but can be important for very low-density mate¬ 
rials close to the BCS-BEC crossover, like recently discov¬ 
ered YPtBi superconductor. For d-wave superconductors 
we calculated both conventional (local) and diffusion- 
induced attenuation rate and show that the latter may 
dominate in strongly anisotropic BSCCO materials at 
moderate temperatures. 
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27r2 

3 Us/ pW 
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where Sk = -|- A^. The same contribution arises 

from the A ( 1 x 2 ) term thus resulting in 


Appendix A: Electron-phonon vertex 

In the main text we have used the tensor-like electron- 
phonon vertex, 

= f {dr)tp(padi3H - {pfVf/ d)5ap^dpUaii^. 

(Al) 

However,careful calculation shows that the actual 
electron-phonon interaction in a co-moving reference 
frame contains additional contribution®. 


(f),^=e(s)(2-Z(e))f3. (A8) 

In general case the vertex Eq.(A8) depends on electron 
energy in a sophisticated way. Meanwhile, in the 
clean limit, when the impurities are negligible and 
e(e) —>■ e, Z{e) —>■ 1, we would have a very simple result 
= ef^. It turns out that this is in fact the correct 
answer once everything is taken into account. 


ttCFR _ ( ttCFR \ I / ttCFR \ 
^el — ph \-^el — ph)^ ' \-^el — ph)‘^'^ 


(A2) 


{Hei-phh= / (dr)V' -C/(r)(a„Wa) V'- 


(A3) 


This contribution describes changes in disorder potential 
strength, for example changes in local concentration of 
impurities. This term is usually negligible. However, 
this is not the case for a d-wave state, where it plays an 
important role. 


2. Full effective vertex 

What we have just left out is the disorder- 

induced contribution to the electron-phonon interaction: 

nodes 


After we average this vertex over the disorder we see that 
1. Naive (incomplete) effective vertex it gives 


Let us first show what would be obtained for an effec¬ 
tive phonon-diffuson vertex ignoring the disorder-related 
contribution. Taking into account only (iL^:0j(^)i we 
have 


(f >df = « E - {pFVF/d)S^p) G^. 

nodes 

(A4) 

For longitudinal phonons, after the summation over 
the nodes is performed, the relevant terms arising from 
the vertex are 


{PadpH - {pFVF/d)6al3) =2 ^Tq -|- A(zf2) 


nodes 


■ 


(A5) 

We will omit Sag and a—,j3— indices for brevity. The 
term for example gives 


(f)dla=“ / (^P) 




{{ek-sY+l^Y 


= e(£) ( 1 - 1 ^3, 


(A6) 


(A7) 


(f)d2 = “ E /(AlO) 

nodes 

— 2ReS'^ = 2^e — e(e)^f3. 


In other words, this vertex is tightly connected to the 
renormalization of electron energy. If energy is not renor¬ 
malized then this vertex is absent. 

Recalling the relation e = Z{e)£{£) we now see that 
the actual effective phonon-diffuson vertex is indeed con¬ 
nected only to the true electron energy variable e: 

(f)^^ = (2e-£)f3, (f)_^2 = (2e-2£>3, (All) 

so that 


(f). = (f).i + (f),2 = "^3, (A12) 

exactly the expression that was used in the main text. 
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